Evidence from d+Au measurements for final-state suppression of high px hadrons in 

Au+Au collisions at RHIC 
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We report measurements of single-particle inclusive spectra and two-particle azimuthal distribu- 
tions of charged hadrons at high transverse momentum (high pr) in minimum bias and central 
d-|-Au collisions at -^snn=200 GeV. The inclusive yield is enhanced in d-|-Au collisions relative to 
binary-scaled p-|-p collisions, while the two-particle azimuthal distributions are very similar to those 
observed in p-l-p collisions. These results demonstrate that the strong suppression of the inclusive 
yield and back-to-back correlations at high pr previously observed in central Au-|-Au collisions are 
due to final-state interactions with the dense medium generated in such collisions. 

PACS numbers: 25.75.-q, 25.75.Dw,25.75.Gz 



Energetic partons propagating through matter are pre- 
dicted to lose energy through induced gluon radiation, 
with the magnitude of the energy loss depending strongly 
on the color charge density |l|. Partonic energy loss is 



potentially a sensitive probe of the matter created in high 
energy heavy-ion collisions, where a quark-gluon plasma 
may form if sufficiently high energy density is achieved. 
The energetic partons originate in the hard scattering of 
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partons from the incoming nuclei. Direct measurement 
of jets resulting from parton fragmentation is difiicult in 
nuclear collisions; nevertheless partonic energy loss can 
be studied using observables such as inclusive spectra 
and two-particle azimuthal distributions of high trans- 
verse momentum (high pt) hadrons. 

Measurements of high pT hadron production in ul- 
trarelativistic interactions of heavy nuclei reveal strong 
suppression of both the single-particle inclusive yield 

IE Q IS and back-to-back pairs (large azimuthal sep- 
aration A(f)) in the most-central, violent collisions, while 
near-side pairs (small Ac^) exhibit jet-like correlations 
that are similar to those in proton-|-proton (p-l-p) col- 
lisions 0. One interpretation of these results is that, 
in the final state following the hard scattering, energetic 
partons traversing the dense medium in the core of the 
collision lose energy, and the observed jets are primar- 
ily those created from partons produced near the surface 
and directed outwards |6| . Alternatively, the suppression 
might result from initial-state effects prior to the hard 
scattering, such as the saturation of gluon densities in the 
incoming nuclei Models incorporating either picture 
are capable of describing central Au-|-Au collision data 
1^. Initial- and final-state effects in Au+Au collisions 
can be separated through studies of deuteron(d)-t-Au col- 
lisions. Theoretical expectations for d-|-Au collisions at 
the Relativistic Heavy Ion Collider (RHIC) are given 
inllliliHIIllIllIllIllI^- Within a perturba- 
tive QCD (pQCD) framework, the expected initial-state 
nuclear effects in d-|-Au collisions are multiple scattering 
prior to a hard collision, which has been used to explain 
the Cronin enhancement of the inclusive yield , and 
shadowing of the parton distribution functions. Nuclear 
effects are expected to increase for more central collisions; 
thus the centrality dependence of observables measured 
in d+Au collisions also will help reveal their origin. 

The STAR Collaboration reports measurements of the 
inclusive invariant pT distribution and two-particle az- 
imuthal distributions at high px for charged hadrons 
[(h+ + h.~)/2, approximated by the summed yields of pri- 
mary TT^, K^, p and p] in minimum bias and central 
d-l-Au collisions at center of mass energy y^snn=200 GeV 
per nucleon pair. Comparison is made to measurements 
at ^snn=200 GeV in the same detector for Au-|-Au and 
p-|-p interactions |E 1^ . The inclusive yield is enhanced 
in d-|-Au collisions relative to binary-scaled p-l-p colli- 
sions, in contrast to the large suppression observed in 
central Au-|-Au interactions. Similar results are reported 
in . The d-|- Au two-particle azimuthal distri- 

butions are very similar to those observed in p-l-p col- 
lisions. These observations are consistent with expecta- 
tions from pQCD models incorporatingboth the Cronin 
enhancement and nuclear shadowing 0, 0, 0, 0, H^ . 
and are inconsistent with calculations that attribute the 
suppression in central Au-I-Au collisions to initial-state 
gluon saturation 0. 
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FIG. 1: Uncorrected charged particle multiplicity distribu- 
tions measured in —3.8 < rj < —2.8 (Au-direction) for d-|-Au 
collisions. Points are for minimum bias (triangles) and periph- 
eral (circles, ZDC-d single neutron) collisions. Both are nor- 
malized to the total number of d-|-Au collisions. Histograms 
are Glauber model calculations. 

STAR is a multi-purpose detector located at 

Brookhaven National Laboratory's RHIC facility. For 
these measurements, the minimum bias trigger required 
at least one beam-rapidity neutron in ZDC-Au, the Zero 
Degree Calorimeter (ZDC) in the Au beam direction, 
which is assigned negative pseudorapidity (ry) . This trig- 
ger accepts 95±3% of the d-l-Au hadronic cross section 
^hadr- Trigger backgrounds were measured using beam 
bunches not in collision. Charged particle momenta were 
measured by the Time Projection Chamber in a 0.5 T 
solenoidal magnetic field. 

After event selection cuts, the data set consists of 10'' 
minimum bias d-|-Au events. Data were analyzed using 
the techniques described in 0, The vertex was re- 
constructed in 93±1% of triggered minimum bias events. 
The spectra were corrected for trigger and vertex-finding 
efficiencies. Contamination of the spectra due to weak 
decay products was corrected based on HI JING . Re- 
sults of an independent analysis, using a different tech- 
nique for vertex reconstruction agree with the re- 
ported spectrum within the relative normalization un- 
certainties at all Pt- 

Centrality tagging of d-f Au collisions is based on the 
raw (uncorrected) charged particle multiplicity within 
—3.8 < rj < —2.8, measured by the Forward Time Projec- 
tion Chamber in the Au beam direction (FTPC-Au |20|). 
The FTPC-Au multiplicity was examined in quadrants 
relative to the orientation of the leading charged hadron 
at mid-rapidity; auto-correlation effects were found to 
be negligible. An independent centrality tag, used as a 
cross-check, requires at least one beam-rapidity (specta- 
tor) neutron in ZDC-d, the ZDC in the deuteron beam 
direction. The cross section for this process in hadronic 
events was measured to be (19.2±1.3)% of cr^f^"^. ZDC-d 
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and FTPC-Au are separated by 8 rapidity units. Figure 
^shows the FTPC-Au multiplicity for minimum bias and 
ZDC-d neutron-tagged events. The latter have a strong 
bias toward low multiplicity. 

The centrality tags were modeled using a Monte Carlo 
Glauber calculation incorporating the Hulthen wave- 
function of the deuteron|2^. In this model the mean 
number of binary collisions (Nbin) is 7.5±0.4 for mini- 
mum bias events and (t^;^^^=2.21±0.09 b. Events with a 



neutron spectator from the deuteron comprise (18±3)% 



of crf^'^dr i'^ model. This event class is biased to- 



ward peripheral coUisions, with (Nbin)=2.9±0.2. The 
FTPC-Au multiplicity distribution was modeled by con- 
voluting the Glauber model distribution of participants 
from the Au nucleus with the charged multiplicity dis- 
tribution measured in 2.5<|77|<3.5 for p-l-p collisions at 
Vs=200 GeVHI. The FTPC-Au acceptance, efficiency 
and backgrounds were taken into account using HIJING 
[21I 1 events in a GEANT model of the detector. Figure 
n shows the measurements for both minimum bias and 
ZDC-d neutron-tagged events, together with the corre- 
sponding Glauber model predictions. The model is val- 
idated by its agreement with both multiplicity distribu- 
tions and with the ZDC-d single neutron cross section 
fraction. High FTPC-Au multiplicity therefore biases 
towards central collisions. Figure ^ shows the cut defin- 
ing the 20% highest multiplicity collisions in the data. 
(Nbin) = 15.0±l.l for the 20% highest multiplicity colli- 
sions in the Glauber model, where the uncertainty in- 
cludes the spread in values obtained with several alter- 
native models. 

Figure 121 shows the invariant inclusive pt distribution 
of (h+ + h^)/2 within |?7|<0.5 for minimum bias and cen- 
tral d-|-Au collisions, together with that for p-f p collisions 
from The error bars are the quadrature sum of sta- 
tistical errors and point-to-point systematic uncertain- 
ties. The normalization uncertainty for d-|-Au collisions 
is 10%. 

Nuclear effects on hadron production in d+Au and 
Au+Au collisions are measured through comparison to 
the p-t-p spectrum using the ratio 



Rab{pt) 



(PN/dprd-q 
Tab d^aPP/dprdri ' 



(1) 



where d^ N / dpTdrj is the differential yield per event in 
the nuclear collision A + B, TAB={^hin) /crf^^i describes 
the nuclear geometry, and d'^a^^ /dprdr] for p-l-p inelastic 
collisions is determined from the measured p-l-p differen- 
tial cross section^. In the absence of nuclear effects 
such as shadowing, the Cronin effect, or gluon satura- 
tion, hard processes are expected to scale with the num- 
ber of binary collisions and Rab{pt) = ^- Figure O shows 
RabIpt) for minimum bias and central d-|-Au collisions. 
The error bars are the quadrature sum of the statistical 
and point-to-point systematic uncertainties. -Rab(j't)>1 
for 2<pt<7 GeV/c. Rab{pt) for central and minimum 
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FIG. 2: Inclusive pr distributions for minimum bias and 
central d+Au collisions, and non-singly diflractive p-l-p colli- 
sions 13 • Hash marks at the top indicate bin boundaries for 
PT >3.8 GeV/c. 
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FIG. 3: Rab{pt) from Eq. 0for minimum bias and central 
d+Au collisions, and central Au+Au collisions Q . The min- 
imum bias d+Au data are displaced 100 MeV/c to the right 
for clarity. The bands show the normalization uncertainties, 
which are highly correlated point-to-point and between the 
two d+Au distributions. 



bias d+Au collisions contain many common uncertain- 
ties, including dependence on the same p+p reference 
spectrum. The ratio of Rab{j>t) for central relative to 
minimum bias collisions, which factors out these com- 
mon uncertainties, is 1.02+0.03 at 4 GeV/c. Rab{pt) 
may be influenced by nuclear shadowing |l3l | and its cen- 
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FIG. 4: (a) Efficiency corrected two-particle azimuthal dis- 
tributions for minimum bias and central d-fAu collisions, and 
for p-|-p collisions 0]. Curves are fits using Eq. 01 with pa- 
rameters given in Table Q (b) Comparison of two-particle 
azimuthal distributions for central d-|-Au collisions to those 
seen in p-l-p and central Au-|- Au collisions |^ . The respective 
pedestals have been subtracted. 



TABLE I: Fit parameters from Eq. |3] Errors are statistical 
only. 





P+P min. bias 


d-f Au min. bias 


d+Au central 


An 


0.081±0.005 


0.073±0.003 


0.067±0.004 


ajv 


0.18±0.01 


0.20±0.01 


0.22±0.02 


Ab 


0.119±0.007 


0.097±0.004 


0.098±0.007 




0.45±0.03 


0.48±0.02 


0.51±0.03 


P 


0.008±0.001 


0.039±0.001 


0.052±0.002 



trality dependence Q. Figure |3| also shows Rab{pt) for 
central Au-|-Au collisions exhibiting large suppression 
in hadron production at high px- 

Figure ^a) shows the two-particle azimuthal distribu- 
tion _D(A0), defined as 



1 1 dN 



trigger 



€ d{A(l)) ' 



(2) 



for minimum bias and central d-|-Au collisions, and for 
p-|-p collisions^. Only particles within |r7|<0.7 are in- 
cluded in the analysis. Ntrigger is the number of particles 
within 4<pT(trig)<6 GeV/c, referred to as trigger parti- 
cles. The distribution results from the correlation of each 
trigger particle with all associated particles in the same 
event having 2 < pT < pT(trig), where e is the tracking 
efficiency of the associated particles. The normalization 
uncertainties are less than 5%. 

The azimuthal distributions in d-|-Au collisions include 



a near-side {A(j) ~ 0) peak similar to that seen in p-l-p and 
Au-|-Au coUisions (6| that is typical of jet production, and 
a back-to-back {A(p ~ tt) peak similar to that seen in p-fp 
and peripheral Au-|-Au collisions that is typical of di- 
jet events. The azimuthal distributions are characterized 
by a fit to the sum of near-side (first term) and back-to- 
back (second term) Gaussian peaks and a constant: 



(b) : D{Acj)) = A 



N- 



2'IT(T]\l 



2TTaB 



-P. (3) 



Fit parameters are given in Table ^ Their systematic 
uncertainties are highly correlated between the data sets, 
and are less than 20% for crjv and less than 10% for all 
other parameters. The only large difference in the az- 
imuthal distributions in p-l-p and d-|-Au collisions is the 
growth of the pedestal P. It increases with increasing 
(Nbin), but is not proportional to (Nbin) as might be ex- 
pected for incoherent production. Both and as ex- 
hibit at most a small increase from p-l-p to central d-|-Au 
collisions. A small growth in as is expected to result 
from initial-state multiple scattering |2J, l25j . The mod- 
est reduction in the correlation strengths and As 
from p+p to central d+Au collisions is similar to that 
seen previously for peripheral Au+Au collisions @. 

Figure 01b) shows the pedestal-subtracted azimuthal 
distributions for p+p and central d+Au collisions. 
The azimuthal distributions are shown also for central 
Au+Au collisions after subtraction of the elliptic flow 
and pedestal contributions . The near-side peak is sim- 
ilar in all three systems, while the back-to-back peak in 
central Au+Au shows a dramatic suppression relative to 
p+p and d+Au. 

The contrast between d+Au and central Au+Au col- 
lisions in Figs. 13 and 0] indicates that the cause of the 
strong high px suppression observed previously is asso- 
ciated with the medium produced in Au+Au but not in 
d+Au collisions. The suppression of the inclusive hadron 
yield at high px in central Au+Au collisions has been 
discussed theoretically in various approaches (see 0] for 
references). Measurements of central Au+Au collisions 
^] are described both by pQCD calculations that incor- 
porate shadowing, the Cronin effect, and partonic energy 
loss in dense matter, and by a calculation extending the 
saturation model to high momentum transfer. However, 
predictions of these models differ significantly for d+Au 
collisions. Due to the Cronin effect, pQCD models pre- 
dict that i?As(pT)>l within 2<pT<6 GeV/c for mini- 
mum bias d+Au collisions, with a peak magnitude of l.f- 
1.5 in the range 2.5<pt<4 GeV/c The enhancement 
is expected to be larger for central collisions ^j. The 
saturation model calculation in predicts RABipT)<^, 
with larger suppression for more central events, achieving 
Rab{pt)^ 0.75 for the 20% most central coUisions. In 
contrast, another saturation model calculation jl5l | gener- 
ates an enhancement in Rab{pt), similar to the Cronin 
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effect, for both d+Au and Au+Au coUisions. Figure El 
shows that Rab{pt) is quahtatively different in d+Au 
and central Au+Au colhsions: in d+Au, Rab{pt) signif- 
icantly exceeds unity. These results are consistent with 
expectations from pQCD calculations but not the satu- 
ration model in . Scattering of the hadronic fragments 
of jets also may contribute to the suppression of the in- 
clusive yield [a, l2y| • 

The azimuthal distributions of back-to-back jets and 
high pt di-hadrons have been observed to broaden in 
fixed-target p+nucleus collisions relative to p+p col- 
lisions, but are not strongly suppressed Slight 
broadening of the back-to-back hadron distribution in 
d+Au collisions at ^snn=200 GeV is also expected from 
pQCD models incorporating the Cronin effect [2^. Pre- 
dictions of the saturation model for the back-to-back 
hadron distributions require further theoretical develop- 
ment, though the rate may be suppressed due to a mono- 
jet contribution Table m shows that the distribution 
of back-to-back high pT hadrons is not substantially mod- 
ified in central d+Au collisions relative to p+p collisions, 
consistent with expectations from pQCD calculations. 

In summary, we have reported the inclusive pt distri- 
butions and two-particle azimuthal distributions of high 
Pt hadrons in minimum bias and central d+Au collisions 
at y^SNN=200 GcV. Similar measurements for Au+Au 
and p+p interactions have revealed a striking suppression 
of both the inclusive hadron yield and the back-to-back 
correlations for central Au+Au collisions. If the suppres- 
sion is the result of initial-state effects, it also should be 
observed in d+Au collisions. No suppression in d+Au 
collisions is observed. Rather, the inclusive yield is en- 
hanced and the two-particle azimuthal distributions ex- 
hibit little change relative to p+p. These results suggest 
that the Cronin effect plays a significant role in d+Au 
collisions for 2 < pT < 7 GeV/c. We conclude that the 
suppression phenomena seen in central Au+Au collisions 
are due to final-state interactions with the dense system 
generated in the collision. 
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